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I . INTRODUCTION 


The  separation  of  a projectile  from  the  gun  tube  releases  the  high 
pressure  propellant  gases  which  can  expand  to  velocities  significantly 
higher  than  the  projectile  velocity.  The  gasdynamic  loadings  generated 
in  the  reverse  flow  can  be  quite  severe  and  possibly  result  in  the  de- 
viation of  the  round  from  its  intended  trajectory.  Since  fin-stabilized 
projectiles  are  designed  to  be  statically  stable  in  forward  flight,  they 
are  unstable  in  the  reverse  flow  near  the  muzzle.  The  present  report 
investigates  the  magnitude  and  duration  of  loadings  experienced  by  a 
fin-stabilized  projectile  in  transit  of  the  muzzle  blast.  A model  is 
developed  to  approximate  these  loads  and  their  effect  upon  the  subse- 
quent trajectory  of  the  round. 

The  most  complete  analysis  of  the  flow  from  the  muzzle  of  a gun 
is  that  of  Oswatitsch^.  In  this  work,  he  examines  the  muzzle  flow 
about  low  and  high  velocity  guns  computing  both  the  one-dimensional, 
unsteady  expansion  which  propagates  into  the  gun  tube  and  the  axially 
symmetric  expansion  of  the  propellant  gases  into  the  atmosphere. 

Based  on  these  calculations,  he  postulates  a quasi-steady  model  of  the 
flow  field  external  to  the  tube,  Figure  1.  Behind  the  unsteady  shock 
layer  advancing  over  the  projectile,  an  underexpanded,  supersonic  jet 
structure  forms  within  the  expanding  propellant  gases.  Oswatitsch 
postulates  that  this  jet  is  quasi-steady  in  the  sense  that  the  core 
properties  are  independent  of  the  unsteady  development  of  the  boundary 
regions  (i.e.,  external  to  the  jet  shocks)  and  vary  only  in  response  to 
changes  in  the  muzzle  exit  conditions. 

This  model  of  Oswatitsch  has  been  experimentally  verified  for 
high  velocity  guns  by  Schmidt  and  Shear^.  They  observe  the  details  of 
the  flow  structure  developing  about  the  muzzle  of  a small  caliber 
rifle  launching  a ball  projectile.  Behind  the  advancing  shock  layer 
an  underexpanded  jet  structure  forms  which  remains  geometrically  in- 
variant throughout  the  period  of  projectile  residence  in  the  muzzle 
blast.  Since  the  muzzle  properties  are  nearly  constant  during  this  time, 
they  conclude  that  the  jet  core  may  be  adequately  represented  by  steady 
analyses . 


1.  K.  Oswatitsch,  "Intermediate  Ballistics,"  Deutsche  Luft  und  Raumfahrt 
FB  64-37 , DVL  Bericht  358 , December  1964.  AD  473249. 

2.  E.  M.  Schmidt,  and  D.  D.  Shear,  "The  Flow  Field  About  the  Muzzle  of 
an  M-16  Rifle,  " BRLR  1692,  U.S.  Army  Ballistic  Research  Laboratories, 
Aberdeen  Proving  Ground , Maryland,  January  1974.  AD  916646L. 


Gretler  uses  the  quasi-steady  model  to  compute  the  loadings  on 
a fin-stabilized  projectile.  Like  Oswatitsch,  Gretler  assumes  the 
highly  underexpanded  jet  core  may  be  approximated  as  a spherical  source. 
The  fin  loadings  are  estimated  using  two-dimensional,  thin  airfoil  theory 
without  correcting  for  flow  inclination,  wing  tip  effects  or  wing-body 
interference.  The  loadings  on  the  projectile  body  are  assumed  negligible. 
These  approximations  completely  eliminate  the  influence  on  wing  and  pro- 
jectile geometry  making  the  model  nearly  universally  applicable.  The 
value  of  loadings  predicted  is  an  upper  bound  on  the  magnitude  of  muzzle 
gas  effects,  thus  providing  a useful  estimate  of  the  possible  influence 
of  muzzle  blast  upon  projectile  trajectory,  i.e.,  dispersion. 

In  the  present  analysis,  the  approach  of  Oswatitsch  and  Gretler  is 
extended  to  include  both  in-bore  gasdynamic  loadings  and  a more  accurate 
model  of  the  muzzle  jet.  The  in-bore  loadings  are  caused  by  passage  of 
the  muzzle  expulsion  over  the  fins.  Figure  2.  This  expansion  occurs 
when  the  propellant  gas  velocity  behind  the  projectile  is  subsonic  prior 
to  launch.  Separation  of  the  obturator  from  the  tube  generates  an  un- 
steady, one-dimensional,  centered  expansion  fan  in  the  propellant  gases, 
expanding  them  to  a sonic  exit  condition.  Through  this  expansion,  the 
flow  velocity  increases  to  values  greater  than  the  projectile  velocity. 

If  the  projectile  is  canted  in  the  tube,  the  flow  over  the  projectile 
generates  transverse  loadings.  The  analysis  of  the  external  muzzle  flow 
field  applies  the  quasi-steady  approach  of  Oswatitsch*;  however,  rather 
than  assuming  a spherical  source  flow,  a more  exact  method  of  character- 
istics computation  of  the  muzzle-jet  flow  is  utilized. 

The  results  of  the  muzzle  flow  calculations  are  used  to  compute 
the  projectile  dynamics  in  this  region.  Of  particular  interest  are 
the  projectile  transverse  velocity,  transverse  angular  velocity,  and 
roll  rate  upon  penetration  of  the  muzzle  blast.  These  quantities  are 
basic  input  into  the  computation  of  the  jump  of  the  round.  The  jump  is 
computed  for  both  symmetric  projectiles  at  angle  of  attack  and  slightly 
asymmetric  projectiles  at  zero  angle  of  attack.  In  the  latter  case, 
the  effect  upon  jump  of  reverse  spin  induced  by  muzzle  gas  loadings  is 
considered. 


II.  MUZZLE  FLOW  ANALYSIS 

Since  fin-stabilized  projectiles  are  generally  launched  with  the 
aid  of  sabots,  it  is  necessary  to  consider  the  impact  of  sabot  con- 
figuration upon  the  muzzle  flow.  The  three  sabot  designs  shown  in 


3.  W.  Gretler,  "Intermediate  Ballistics  Investigations  of  Wing 
Stabilised  Projectiles,"  German  Air  and  Space  Research  Report 
67-92,  FSTC-HT- 23-22-69-7 2,  1967. 
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Figure  3 are  arranged  in  order  of  increasing  separation  between  the 
sabot  and  fin  assemblies.  The  cup  or  pusher  sabot  completely  surrounds 
the  fins.  Obviously,  the  muzzle  gasdynamics  of  this  design  would 
resemble  those  of  a full-bore,  ball  projectile.  The  data  of  Schmidt 
and  Shear  show  that  a ball  projectile  is  not  immersed  in  the  quasi- 
steady core  flow,  but  rather  interacts  strongly  with  the  developing 
shock  layer.  This  flow  is  dominated  by  unsteady  effects  and  is  diffi- 
cult to  treat  analytically.  However,  the  cup  sabot  is  generally  used 
to  launch  spin-stabilized  projectiles,  which  are  not  within  the  scope 
of  the  present  survey. 

Fin-stabilized  projectiles  employ  ring  or  puller  sabots  to  facili- 
tate launch.  With  these  designs,  the  projectile  fins  extend  behind  the 
sabot.  Sufficient  extension  eliminates  aerodynamic  interference  be- 
tween the  fins  and  sabot,  permitting  direct  exposure  of  these  strong 
lifting  surfaces  to  the  reverse  muzzle  flow.  For  a flechette  round, 

Glauz  shows  that  the  most  significant  transverse  loadings  on  the  pro- 
jectile are  generated  at  the  fin  surfaces.  He  further  notes  that  a 
first  approximation  of  the  transverse  loadings  may  be  computed  by  con- 
sidering solely  the  fin  forces. 

To  illustrate  the  development  of  the  flow  about  the  fins  during 
shot  ejection,  consider  the  diagram  in  Figure  4.  Since  the  obturator 
is  located  forward  on  the  projectile,  the  propellant  gas  is  released 
while  the  fins  are  still  in-bore.  For  subsonic  propellant  gas  velocities, 
one-dimensional  expansion  waves  propagate  back  up  the  tube,  accelerating 
the  gases  to  a sonic  velocity  at  the  muzzle.  The  initial  projectile 
and  propellant  gas  velocities  are  equal;  thus,  passage  of  the  expansion 
waves  over  the  fins  results  in  a reverse  flow  relative  to  these  surfaces, 
increasing  from  zero  to  a maximum  value  at  the  muzzle.  External  to  the 
muzzle,  the  underexpanded  propellant  gas  jet  and  free  air  blast  develop. 
TJie  region  between  the  muzzle  and  the  Mach  disc  forms  the  quasi-steady, 
supersonic  core  of  the  jet.  Since  the  Mach  disc  is  attached  to  or 
precedes  the  obturator  on  the  sabot  base,  the  fins  are  exposed  to  the 
jet  core.  The  propellant  gas  velocity  increases  rapidly  from  the  muzzle, 
causing  the  relative  velocity  over  the  fins  to  vary  from  subsonic  through 
supersonic  values.  This  direct  exposure  of  projectile  fins  to  the  super- 
sonic jet  core  has  been  experimentally  observed  . A spark  shadowgraph 
of  the  muzzle  flow  about  a 5.77mm  smoothbore  gun  firing  a saboted. 


4.  W.  D.  Glauz 3 "Estimation  of  Forces  on  a Flechette  Resulting  from 
a Shock  Wave , " Midoest  Research  Institute , Kansas  City , Missouri , 
Final  Report , Project  No.  S451-E,  May  1971.  AD  724178. 

5.  E.  Mo  Schmidt  and  D.  D.  Shear , "Discard  of  XM-645  Sabot  in  Muzzle 
Blast>"  BRL  Memorandum  to  be  published  in  1975. 
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flechette  round  is  shown  in  Figure  5.  The  rear  of  the  projectile  is 
clearly  outlined.  From  the  observed  position  of  the  Mach  disc,  it  is 
apparent  that  the  fins  are  well  within  the  supersonic  core  of  the  muzzle 
jet. 


As  the  outward  propagation  of  the  muzzle  gases  slows  due  to  radial 
expansion,  the  projectile  penetrates  the  shock  layer  and  enters  free 
flight.  Compared  with  the  intensity  of  loadings  seen  near  the  muzzle, 
the  forces  exerted  cn  the  projectile  in  the  far  field  of  the  jet  and 
in  the  shock  layer  are  negligible.  Oswatitsch1  notes  that  within  a 
few  calibers  of  the  muzzle  the  thrust  exerted  on  a projectile  drops  by 
many  orders  of  magnitude.  Schmidt  and  Shear^  use  a semi-empirical 
approach  to  determine  property  values  in  the  shock  layer,  i„e.,  the 
region  between  the  Mach  disc  and  air  blast.  Their  analysis  predicts 
pressure  levels  in  the  shock  layer  of  3-4  atmospheres  or  orders  of 
magnitude  less  than  the  pressures  near  the  muzzle  which  reach  several 
hundred  atmospheres. 

In  the  remainder  of  this  section,  techniques  will  be  developed  for 
the  calculation  of  a first  order  estimation  of  the  gasdynamic  loadings 
experienced  by  a fin-stabilized  projectile  during  launch.  It  will  be 
assumed  that  the  ring  or  puller  sabot  design  is  being  used  and  that 
there  is  no  interaction  between  the  projectile  and  the  sabot.  Following 
Glauz4,  we  consider  only  the  contribution  of  fin  surfaces  to  the 
generation  of  transverse  loads.  Both  the  in-bore  and  external  flow 
fields  will  be  considered.  The  exterior  flow  will  be  modelled  using 
the  quasi-steady  approximation  of  Oswatitsch.  No  attempt  will  be  made 
to  include  the  loadings  generated  in  transit  of  the  shock  layer.  Two 
subsections  will  be  presented.  The  first  deals  with  development  of 
the  free  jet  model  and  fin  loadings.  The  second  presents  a similarity 
analysis  of  the  in-bore  expansion. 

A.  Flow  Exterior  to  the  Muzzle 

The  analysis  in  this  section  estimates  the  fin  loadings  during 
transit  of  the  propellant  gas  jet.  No  attempt  is  made  to  compute  the 
effect  of  passage  of  the  fins  through  the  shock  layer.  The  analyses 
of  Oswatitsch  and  Gretler  use  a spherical  source  model  of  the  flow 
within  the  underexpanded,  propellant  gas  jet.  While  such  an  approxi- 
mation is  adequate  ir.  the  far  field,  it  does  not  accurately  represent 
the  flow  properties  in  the  vicinity  of  the  muzzle.  In  Figure  6a,  the 
centerline  Mach  number  distribution  predicted  by  Oswatitsch  is  compared 
with  a more  exact  calculation  obtained  by  applying  the  method  of 
characteristics6.  The  source  flow  model  shows  a very  rapid  expansion 


6.  A.  B.  Vickj  E.  H.  Andrews t J.  S.  Dennardj  and  C.  B.  CraidoHj 

"Comparison  of  Experimental  Free-Jet  Boundaries  with  Theoretical 
Results  Obtained  with  the  Method  of  Characteristics , " NASA  Technical 
Note  D~2Z273  vune  1964.  (National  Technical  Information  Service 
N64-2Z0Z2) 


of  the  flow  near  the  muzzle  due  to  the  constant  radial  divergence 
imposed  by  the  assumption  of  spherical  symmetry.  In  contrast,  the 
method  of  characteristics  computation  shows  a more  gradual  expansion 
of  the  flow  through  a two-dimensional,  Prandtl-Meyer  expansion  at  the 
muzzle.  Since  this  is  a more  realistic  model  of  the  flow,  the  results 
of  the  method  of  characteristics  computation  are  adopted  in  this  analysis. 


Owen  and  Thornhill7  show  that  the  flow  within  the  bounding  shocks 
of  an  underexpanded  free  jet  is  universal  in  nature;  i.e.,  the  flow 

i . . i r»  . • _ _ _ f A.  T_  » • rx  rtm  rt  -*%  "1  Z' 


parameters,  P/Pe>  P/Pe>  T/Te>  V/Qe>  are  functions  of  their  geometric 


locations  within  the  jet,  the  exit  Mach  number,  and  the  ratio  of  specific 


AVI  V.av.XVJU  -t.  ^ J ^ ~ f ' _ . . 

heats,  but  they  are  independent  of  the  jet  pressure  ratio,  Pg/P^-  Since 


the  ratio  of  specific  neats  of  propellant  gases  has  a limited  range 
(1.20  <y<  1.30),  a representative  value,  y=1.2S,  is  selected  to  be  valid 
for  all  type  propellants.  The  Mach  number  of  the  propellant  gases  be- 
hind the  projectile  may  vary  from  subsonic  through  supersonic  values 
depending  on  the  weapon  system  considered.  In  the  subsonic  case,  an 
unsteady  expansion  fan  propagates  into  the  gun  tube  bringing  the  muzzle 
exit  velocity  to  a sonic  value.  Thus,  only  sonic  or  supersonic  exit 
conditions  need  be  considered.  Centerline  Mach  number  distributions 
predicted  by  the  method  of  characteristics  for  supersonic  exit  conditions 
are  compared  with  the  sonic  exit  distribution  in  Figure  6B.  The  flow 
expansion  starts  immediately  downstream  of  the  muzzle  for  the  sonic 
exit  condition;  however,  with  supersonic  exit  velocities,  the  expansion 
waves  from  the  muzzle  lip  are  inclined  in  the  downstream  direction.  The 
centerline  flow  does  not  sense  the  change  in  surroundings  until  these 
waves  arrive  on  axis.  Once  started,  the  centerline  expansion  occurs 
more  rapidly  in  the  case  of  the  supersonic  exit  conditions. 


The  method  of  characteristics  calculation  of  the  jet  presumes  steady 
To  account  for  temporal  variation,  the  quasi-steady  approximation 

_ . _ . . • • « . r -i  4-  U 


postulates  that  the  flow  field  properties  vary  instantaneously  with 
changes  in  the  exit  conditions.  This  approximation  implies  that  the 
signal  propagation  velocity  through  the  muzzle  jet  is  infinite.  However, 
rather  than  account  for  temporal  variations,  Gretler’  assumes  the  exit 
properties  remain  constant  (at  the  values  reached  after  1-D,  unsteady 
expansion  into  the  tube)  throughout  his  calculation.  Schmidt  and  Shear 
show  the  muzzle  pressure  to  be  nearly  invariant  during  projectile  residence 
in  the  muzzle  jet  of  high  velocity,  long  tube  ~uns . . Additionally,  in 
moving  through  the  region  of  primary  interest  (within  a few  calibers  of 
the  muzzle  where  loadings  are  maximum)  the  projectile  will  experience 
spatial  property  gradients  orders  o^  magnitude  greater  than  temporal 
variations.  For  these  reasons,  the  subsequent  analysis  of  fin  loadings 


7.  P.  Owen  and  C.  Thornhill > " The  Flow  in  an  Axially  Symmetric  Super- 
sonic Jet  from  a Nearly  Sonic  Orifice  into  a Vacuum , " RARDE , 
Report  30/48,  1948.  AD  57261. 
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in  the  muzzle  jet  will  assume  constant  exit  properties. 


The  coordinate  systems  for  these  calculations  are  shown  in  Figure  7. 
They  were  selected  in  accordance  with  established  coordinates  used  in 
the  analysis8'9  of  projectile  trajectories.  The  nonrolling  coordinates 
are  used  in  the  solution  of  the  projectile  angular  motion,  while  the 
earth  fixed  coordinates  are  the  system  in  which  projectile  displacement, 
i.e„,  jump,  is  determined.  In  the  vicinity  of  the  muzzle,  it  will  be 
assumed  that  the  projectile  velocity  vector  lies  along  the  Xg  axis, 

which  is  also  the  axis  of  symmetry  of  the  muzzle  jet  flow. 

Since  the  projectile  attitude  at  exit  is  somewhat  arbitrary,  it 
may  be  assumed,  without  loss  of  generality  that  the  projectile  is 
launched  with  an  initial  angle  of  attack  but  zero  sideslip.  The  initial 
transverse  linear  and  angular  velocities  of  the  projectile  imparted  by 
in-bore  mechanical  loadings  are  unknown;  however,  since  the  time  scales 
under  consideration  are  quite  small,  it  will  be  assumed  that  these 
velocities  do  not  result  in  significant  variations  in  projectile 
attitude  in  the  muzzle  flow  (measurements  show  only  a few  mil  variation 
for  projectiles  with  first  maximum  yaws  of  10°  or  less)  . Since  the 
values  of  transverse  velocities  at  the  muzzle  do  not  enter  into  the 
calculation  of  fin  loadings  in  the  propellant  gas  jet,  it  is  convenient 
to  assume  them  to  be  zero.  However,  in  a calculation  of  all  transverse 
loadings,  the  .in-bore  and  muzzle  blast  effects  could  be  summed  vector- 
ially . 

The  lift  on  fin  surfaces  due  to  differential  fin  cant  generates 
a rolling  moment,  while  the  lift  due  to  projectile  angle  of  attack  pro- 
duces a tran  verse  force  and  overturning  moment.  The  lift  on  a single 
fin  due  to  muzzle  jet  flow  may  be  expressed  as: 

* - C£  1 “Vr  A (1) 

where:  C£  - Lift  Coefficient 

p - Local  (Jet)  Density 
A - Fin  Planform  Area 

and  V is  the  relative  velocity  between  the  projectile,  V , 
r P 


8.  C.  H.  Murphy,  "Free  Flight  Motion  of  Symmetric  Missiles,  " BRLR  1216 
U.S.  Army  Ballistic  Research  Laboratories,  Aberdeen  Proving  Ground, 
Maryland,  July  1963.  AD  442757. 

9.  C.  H.  Murphy  and  J,  W.  Bradley,  "Jump  Due  to  Aerodynamic  Asymmetry 
of  a Missile  with  Varying  Roll  Rate,"  BRLR  1077,  U.S.  Army  Ballistic 
Research  Laboratories,  Aberdeen  Proving  Ground,  Maryland,  May  1959. 
AD  219312. 
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and  the  local  jet  flow,  V: 


V = V - V . 
r p 


For  small  angles  of  attack,  the  lift  on  the  fin  will  be  assumed 
to  be  a linear  function  of  the  angle  of  attack  of  the  fin: 


CjC  ~ CjC  + CjC  a • 


In  the  absence  of  camber,  the  value  of  C£  is  zero.  will  be 

0 a 


evaluated  using  two-dimensional,  linearized  airfoil  theory.  This 
approach  neglects  the  effects  of  wing  tips,  wing-body  interference, 
and  the  presence  of  a stem  shock  on  the  projectile.  Following  the 
approach  of  Gretler'  , the  Prandtl-Glauert  Rule  will  be  used  for  sub- 
sonic flow. 


-=■  For  0 < M < 0.7, 
2 r — * 


an  .2  Ackeret  Airfoil  Theory  for  supersonic  flow 


For  M > 1.1, 
r — ’ 


where 


The  value  of  approaches  infinity  in  both  equations  (4)  and  (5) 
a 

as  the  relative  Mach  number,  M , approaches  one.  The  specific  bounds 


on  these  relations  are  dependent  upon  the  angle  of  attack  and  thinness 
of  the  airfoil.  For  thin  airfoils  and  small  angles  of  attack,  the 
relations  may  be  used  quite  close  to  the  sonic  condition.  To  obtain  an 
epproximate  coefficient  in  the  transonic  region,  a constant  value^  of 
C j,  will  be  assumed: 
a 

Cj,  =8.8  For  0.7  < M <1.1  (7) 

a r 


The  resulting  behavior  of  the  lift  coefficient  with  increasing  Mach 
number  is  shown  in  Figure  8. 


Using  this  linear  coefficient,  the  lift  on  a single  fin  is: 


* ■ c£  a T pvr2  A 

a 


where:  p*,  p*,  T*  refer  to  critical  (or  sonic)  values  of  the  flow 

properties,  and:  y 


M * = 
r 


/yr'F' 


(10) 


Introducing  the  critical  momentum  : 


* 4.  n*  r*  B (y+1)  p*  , 


p*  + p*  c 

a non-dimensional  lift  function  is  defined  as 


(ID 


L = 


(p*  + p*  c*  ) AS 


(Y+l)  p*  Aa 


= C, 


Y 

2(Y+lT 


£t  m 
p*  r 


(12) 


Since  the  flow  under  consideration  is  isentropic,  the  expression 
on  the  right  hand  side  of  Equation  (12)  may  be  reduced  to  a function 
of  the  jet  Mach  number  and  a non-dimensional  projectile  velocity, 

V /c  , where  c^  is  the  speed  of  sound  of  the  propellant  gas  prior  to 

shot  ejection.  Since  the  projectile  velocity  is  taken  as  constant 
through  the  muzzle  gases,  the  parameter  Vp/c1  is  also  the  propellant 

gas  Mach  number  prior  to  shoe  ejection.  Using  the  isentropic  relations 


(—  + Hi  M2)"y-1 

‘‘Y+l  Y+l 


= M 


(13) 

(14) 


r c’ 


— ar 
C* 


(JL  + Hi  M2) 
'•y+l  Y+1 


(15) 


and 


V 


M 


- Vp  Vp 

c r 


v C. 

= . _E  1 £l 

C1  c*  c 


(16) 


For  subsonic  propellant  gas  velocities  prior  to  shot  ejection 
(V  /c,  < 1.0),  one-dimensional,  unsteady,  isentropic  flow  theory 
v p 1 
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predicts: 


I 


[ 


1 - t 2 n + Y_1 
^ l y+T  ~T~ 


v 

a 


-l 


(17) 


Also: 


V 

A 


2y 

Y-l 


(18) 


For  sonic  propellant  gas  velocities  prior  to  shot  ejection, 

Cj  - c*,  and  for  supersonic  velocities  Equation  (15)  may  be  used  to  com- 
pute c^/c*  by  substituting  V /c^  for  M. 


As  previously  stated,  substitution  of  Equations  (13)  - (17)  into 
Equation  (12)  results  in  an  expression  for  the  non-dimensional  lift 
coefficient  which  is  a function  of  M and  V /c,  only.  Since  the  center- 

p 1 

line  Mach  number  distribution  is  given  in  Figure  6,  the  lift  function, 

L,  may  be  evaluated  along  the  jet  axis  for  various  values  of  V /c, , 

. — P 1 

Figure  9.  The  variation  of  L at  X /D  = 0 with  V /c„  reflects  the 

e p 1 

effect  of  the  in-bore  expansion.  For  values  of  V /c  less  than  one,  the 

expansion  accelerates  the  propellant  gases  to  velocities  greater  than 
the  projectile  velocity;  thus,  at  the  muzzle,  X /D*0,  th''  projectile 


experiences  a velocity  differential  which  increases  with  the  strength 
of  the  expansion  (i.e.,  with  decreasing  V^/Cj) . For  values  of 

Vp/c^  >_  1.0,  the  relative  velocity  between  the  projectile  and  pro- 
pellant gases  remains  zero  until  they  experience  the  two-dimensional, 
external  expansion. 


The  initial  increase  and  subsequent  decay  of  the  lift  function 
reflects  the  property  variations  in  the  free  jet  flow  field  and  their 
relation  to  the  parameters  in  Equation  (12) : 


L = C 


2(y+1)  P 


P M*2 


The  centerline  property  distribution  of  the  muzzle  jet  considered  is 
shown  in  Figure  10,  Consider  the  case  of  sonic  propellant  gas  velocity, 
Vp/c^  = 1.0.  At  the  muzzle,  the  relative  velocity  is  zero.  Through 

the  jet  expansion  field,  the  propellant  gas  velocity,  V/c*,  increases. 
The  projectile  velocity  remains  constant;  thus,  M^*  increases  with  the 

jet  velocity  causing  L to  increase.  While  the  propellant  gas  velocity 
increases,  the  static  thermodynamic  properties  of  the  jet  flow  rapidly 
decrease,  V/c*  increases  less  rapidly  after  two  calibers,  while  the 
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density  continues  to  decrease,  resulting  in  a rapid  decay  of  the  lift 
function.  Thus,  in  the  first  two  calibers  of  travel  through  the  muzzle 
jet,  the  fins  experience  the  most  significant  part  of  the  gasdynamic 
loadings . 

The  lift  in  itself  is  not  of  primary  interest;  rather,  the  effect 
of  this  lift  upon  the  projectile  dynamics  is  being  sought.  To  determine 
this  effect,  the  impulsive  momentum  transferred  to  the  projectile  must 
be  determined.  For  the  single  fin  considered  thus  far,  Newton's  law 
may  be  written: 


dP 

Ht 


(19) 


where  P;  is  the  momentum  imparted  to  the  projectile  by  a single  fin*. 
The  lift  may  be  integrated  to  obtain  the  change  in  momentum  through 
the  muzzle  blast: 


where 


Pj  « / £ dt 

tl 

= 0 (by  assumption) 
t1  = Time  fins  pass  muzzle 


(20) 


t = Time  at  which  momentum  pulse  is  to 
be  evaluated 

Since  the  lift  is  evaluated  as  a function  of  position 
along  the  jet  axis,  it  is  necessary  to  transform  the  time  variable 
into  projectile  displacement: 


where 


Thus  : 


dt  = — dX 
P 


X 


X 

_e 

D 

D 


P = Tf—  f 
V o 

P 


(21) 


jC  dX 


(22) 


Substituting  the  non-dimensional  lift  function.  Equation  (12),  into 
Equation  (22) : 


p = f (y+i)  p*Aaf  L 


dX 


(23) 


*White  vector  notation  is  not  used  jC  and  P are  vector  quantities . 
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The  integral  on  the  right  hand  side  is  defined  as  the  dimensionless 
momentum  function: 

P = /*  L dX  (24) 

o 

The  evaluation  of  this  integral  follows  directly  from  the  lift  function, 
Figure  9;  values  at  various  jet.  stations  are  given  in  Figure  11. 


The  momentum  function  increases  rapidly  over  the  first  caliber; 
however,  decay  of  the  lift  function  greatly  reduces  the  moment  transfer 
during  the  remainder  of  projectile  travel.  A nearly  asymptotic  value 
of  momentum  transfer  is  obtained  after  only  five  calibers  of  travel^ 

The  limit  of  the  asymptote  is  the  total  momentum  imparted  by  the  je,., 

P , and  is  presented  as  a function  of  V /Cj  in  Figure  12,  For 

V /c  < l.O.  P varies  only  slightly;  however,  for  V /c  > 1.0  P 
pi—  P 1 

rapidly  decays.  The  formulation  of  the  aerodynamics  is  based  on  linear, 
two-dimensional  airfoil  theory  with  no  consideration  given  to  wing  or 
body  geometry;  thus,  while  admittedly  approximate,  the  results. 

Figures  9,  11,  and  12,  are  applicable  to  any  finned  projectile  tra- 
versing the  muzzle  jet  flow  field. 


The  lift  and  momentum  functions  presented  are  for  a single  fin 
of  a projectile.  The  combined  effect  of  multiple  fins  is  evaluated 
by  vector  summation  of  the  impulses  on  the  individual  fins.  Two  con- 
figurations are  of  immediate  interest,  a projectile  with  differentially 
canted  fins  and  a projectile  at  angle  of  attack.  Figure  13  Since  lin- 
ear theory  is  used,  the  two  effects  are  separable  and  will  be  con- 
sidered individually.  Differential  fin  cant  is  used  to  produce  a 
rolling  moment  which  induces  spin  in  free  flight;  however,  in  the 
muzzle  flow  field,  the  direction  of  the  moment  is  reversed.  For  pro- 
jectiles which  are  unspun  at  launch,  the  muzzle  jet  moment  will  produce 
a reverse  spin,  which  can  have  a significant  impact  on  jump  due  to 
asymmetries . 

The  present  model  is  readily  applied  to  the  computation  of  the 
change  in  spin  (roll-rate)  through  the  muzzle  jet.  Assuming  all  fins 
to  be  canted  at  the  same  orientation  to  the  projectile  axis,  the  total 
moment  acting  on  the  projectile  is: 

M,  = - n r £ 

$ o 


= I 24  (25) 

x dt2 

n - Number  of  fins 

<j,  - Roll  angle  (positive  in  the  sense  of  right 


where 


hand  spin) 


r - Distance  from  projectile  axis  to  fin  center  of 


pressure  (taken  as  centroid  of  area). 


Equation  (25)  may  be  integrated  to  obtain  the  change  in  spin  through 
the  muzzle  jet: 

nr  t 


* - J = - — ° / 0 £ dt,  (26) 

To  1 I * 

X tx 


where 


= roll  rate  at  muzzle. 


The  roll  rate  is  transformed  to  the  angular  roll  per  caliber  of 
travel : 


= _R  4*  = _R  / 

dt  l ds  £ 


where 


£ = projectile  shaft  diameter 
X 


Substituting  Eauation  (27)  and  using  Equations  (20)  - (24)  in  Equation 
(26), 


/ / 


- *1 


O D£  , , . * A i n 

j — -J  (y+1)  p*  a A Pc 


x V 


This  expression  permits  the  evaluation  of  the  spin  imparted  to  any  pro- 
jectile in  transit  of  the  muzzle  blast  if  the  projectile  physical 
properties  and  gun  launch  conditions  are  known. 


The  second  configuration  under  consideration  is  a projectile 
traversing  the  muzzle  jet  at  anjle  of  attack,  Figure  13B.  The  transverse 
force  acting  on  a projectile  with  n fins  can  be  equated  to  the  force  on 
a two-finned  projectile  having  the  fins  perpendicular  to  the  plane  of 
the  diagram  and  at  the  angle  of  attack  of  the  projectile.  In  Appendix 
A,  the  combined  plan  form  area  of  the  equivalent  fins  is  shown  to  be 
n/2  times  the  area  of  a single,  original  fin  (for  n>3) . The  lirt  or  transverse 
force  transmitted  to  the  projectile  by  the  fin  assembly  nay  be  expressed 
using  Equation  (12), 


L = (y+1)  p*  -y  A a L 


■ . '• 


, 


»r.  ^sd, 


p at2 

where  A - Area  of  a single  fin 

M - Mass  of  the  projectile 
P 

The  moment  acting  about  the  projectile  center  of  gravity  due  to  this 
lift  force  is: 


M = L (A  cos  a) 


s L A 


(For  small  5) 


The  transverse  linear  and  angular  velocities  imparted  to  the  projectile 
in  the  muzzle  blast  may  be  obtained  by  integrating  Equations  (29)  and 
(30).  The  transverse  linear  velocity  is: 

f - f0  (31) 

d D MV 

P p P P 

and  the  transverse  angular  velocity  is: 

d'  - d'  = (Y+l)  P*  \ A 5 * p0  t32) 

y p 

Again,  the  subscript  one  refers  to  muzzle  properties. 

The  relations  developed  in  Equations  (28) , (31) , and  (32)  permit  the 
evaluation  of  projectile  dynamics  during  transit  of  the  muzzle  blast.  In 
the  next  subsection,  the  effect  of  the  in-bore  expansion  is  considered. 

The  initial,  free  flight  dynamics  of  the  round  will  be  related  to  its 
subsequent  trajectory  in  Section  III  and  IV.  Prior  to  leaving  this  sub- 
section,  the  results  of  the  present  model  of  muzzle  blast  loadings  will  be 
compared  against  previous  approaches . 

The  major  modification  incorporated  is  the  use  of  a method  of  3 
characteristics  calculation  of  the  quasi-steady  core  flow.  Gretler  makes 
direct  use  of  the  spherical  source  model  of  Oswatitsch1.  more  rapid 

flow  expansion  predicted  by  the  source  model  is  reflected  in  the  resul- 
ting lift  and  momentum  functions.  Figure  14.  The  lift  function  computed 
by  Gretler  reaches  a maximum  and  decays  significantly  sooner  than  the 
present  model.  Since  the  momentum  function  is  the  area  under  the  lift 
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curve,  the  momentum  function  predicted  by  Gretler  is  only  half  of  that 
obtained  with  the  present  approach.  These  plots  demonstrate  the  im- 
portance of  accurately  depicting  the  flow  near  the  muzzle  where  maximum 
loadings  occur.  For  more  complex  muzzle  phenomena,  e.g.,  ball  pro- 
jectiles, it  may  be  necessary  to  consider  fully  three-dimensional,  time 
dependent  models  of  the  muzzle  exit  flow. 


B. 


Flow  in  the  Gun  Bore  Subsequent  to  Separation 


1 


Subsonic  velocity  in  the  propellant  gases  driving  the  projectile 
allows  an  expansion  fan  to  propagate  into  the  gun  barrel  subsequent  to 
separation.  The  expansion  accelerates  the  gases  to  a sonic  velocity  at 
the  muzzle,  producing  flow  velocities  within  the  fan  which  are  greater 
than  the  projectile  velocity.  The  relative  flow  over  the  fins  imparts 
momentum  to  the  projectile.  If  the  round  is  canted  with  respect  to  the 

gun  bore  (but  not  touching  the  walls)  transverse  momentum  will  be  gen- 
erated. 


The  analysis  of  the  in-bore  flow  is  simplified  by  assuming  the 
expansion  over  the  projectile  is  one-dimensional.  Neglecting  two- 
dimensional  effects  near  the  muzzle  and  area  change  due  to  the  projectile 
configuration,  the  flow  may  be  represented  as  a centered  expansion. 

Landau  and  Lifshitz  u demonstrate  that  the  centered  expansion  is  a 
similarity  flow;  i.e.,  the  flow  pattern  does  not  change  with  time  if 
the  following  transformed  special  variable  is  used: 

x/t  = V - c . (33) 

With  this  independent  variable,  the  equations  of  motion  may  be  in- 
tegrated. 


V = 


V 

P 


Substitution  of  the  isentropic 

P = 


relation. 

JL 


-vY-1 

Sc  ' 

L 


yields 


V = 


V 

P 


d(c  - c) 

+ I 

Y - 1 


(34) 


(35) 


(36) 


10 o L,  D.  Landau  and  E.  M.  Lifshitz , Fluid  Mechanics , Tevgccmon  Vvess , 
Addison-Wesley  Publishing  Company , Inc. , Reading,  Mass.,  1959  , 
pp  353-359, 
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(37) 


Using  Equation  (33)  to  eliminate  the  speed  of  sound  in  Equation  (36) 

V = [ 2 | + (Y-l)  Vp  + 2 Cj]  / (y+1)  , 

or,  eliminating  the  flow  velocity, 

c = [-(Y-l)  | + (Y-l)  Vp  + 2 Cl]/(Y+1).  (38) 

These  equations  describe  the  flow  seen  by  the  fins  as  they  exit  the 
gun  tube.  The  trajectory  of  the  center  of  pressure  (area  centroid)  of 
the  fins  is 

X = -d  + V t , (39) 

e P 

where  d = distance  from  obturator  to  c.p.  of  fins. 

Rearranging  Equation  (39) : 


XXX 

t'  = vp  (/)/d  ♦ /)  . (40) 

permits  substitution  into  Equations  (37)  and  (38)  to  obtain  expressions 
for  V/c1  and  c/c1  as  functions  of  fin  location,  fin  setback,  and  pro- 

pellant  Mach  number  prior  to  shot  ejection,  Vp/Cj.  Knowing  the  flow 

velocity  and  speed  of  sound,  it  is  now  possible  to  compute  fin  loadings. 

The  lift  function,  Equation  (12),  is  again  applied: 

£.  M*2 

L - C£  Tfw)  p*  r 
a 


I 


f 

I 


T 


Through  an  unstean>  flow  field,  the  local  critical  conditions  are  not 
necessarily  constant;  thus,  the  starred  quantities  in  Equation  (12*  ) are 
evaluated  at  the  muzzle  and  held  constant.  The  relative  Mach  number,  M*> 

is  defined  in  Equations  (14)  and  (16): 


= (M 


V c,  * 

p 1 C*'  £ 
C^  C*  C ^ C* 
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l'he  local  Mach  number  is  simply  the  ratio  of  Equations  (37)  and  (38)  : 


(2  (VC1}  (Xe/d)/(1  + Xe/d)1  + 1 (Y_1)  iVv/ci]  + 21 

M = *■  • “ (41) 

(-(Y-l)  (Vp/C;L)  (Xe/d)/(l  + Xe/d)]  + [ (Y-l)  (V  /Cj)  + 2]  . 

c /c*  is  defined  by  Equation  (17)  and  c*/c  is  derived  from  Equation 
(38): 

c*  (Y-l) (Vp/Cp)  +2 

c = (-(Y-l)  (Vp/cp  (Xe/d)/a  + Xe/d)]  + [ (y-ID  (V  /c1)  + 2]  . (42) 

The  value  of  p/p*  in  the  lift  function  is  obtained  from  the  isentropic 
relation : 

. 1— 

P.  _ ,£%  Y"1 

P*  " lc  J . (43) 

The  variation  of  lift  coefficient,  C , with  Mach  number  is  again 

a 

taken  from  Figure  8.  The  portion  of  the  plot  used  for  the  in-bore 
calculations  is  limited  to  Mach  numbers  between  zero  and  one.  In  this 
range,  there  is  reason  to  question  the  validity  of  linear  theory  to 
model  the  flow  over  the  fins  within  the  confined  geometry  of  the  gun 
tube.  Thus,  the  results  of  this  portion  of  the  report  should  be  viewed 
as  an  initial  attempt  to  define  the  magnitude,  duration,  and  relative 
importance  of  in-bore,  gasdynamic  loadings  during  launch. 


The  momentum  impulse  function  is  determined  by  integrating  the  lift 
function  over  the  fin  trajectory: 

Pb  = V (y+d  p*  M^°  1 dX 

p -X 

r w 


The  non-dimensional  momentum  function  is  defined  as: 


P,  - f°  L dX 
b - X 


The  lower  limit  of  integration  is  the  axial  location  at  which  the  first 
wave  of  the  expansion  fan  intersects  the  fin,  which  is  readily  deter- 
mined to  be : y 

x =-1  (1  + -E)_1  . 

w D v c/ 
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Permitting  the  integral  in  Equation  (45)  to  be  expressed  as 

X 


■»/ 


L d(-|) 


- ci‘V'irl 


The  integral  in  Equation  (46)  is  independent  of  fin  setback  distance,  d, 
varying  only  as  a function  of  the  parameter  V /c^;  thus,  the  momentum 

impulse  function  has  a linear  dependence  on  d/D.  Evaluation  of  for 

d/D  - 1 and  various  V^/c^  permits  the  calculation  of  P^  for  any  d/D  by 

linear  extrapolation,  i.e.: 


= n P, 


(d/D  «>i  ) 


’(d/D  = 1) 


The  in-bore  momentum  impulse  function  variation  with  V.  /c^  is 

compared  to  the  muzzle  jet  function  in  Figure  15.  Not  surprisingly,  the 
in  bore  momentum  function  decreases  rapidly  with  increasing  propellant 
gas  Mach  number,  V /c^.  At  V^/c^  *»  1.0,  no  waves  propagate  into  the  gun 

tube,  and  the  transverse  loadings  can  not  occur.  High  velocity  guns 
generally  have  propellant  gas  Mach  numbers  greater  than  0.5;  thus,  in- 
bore fin  loadings  may  be  neglected.  However,  low  velocity  guns  such  as 
mortars  can  have  very  low  pre- launch,  propellant  gas  Mach  numbers,  and 
in-bore  loadings  will  be  significant. 


III.  DEFLECTION  OF  SYMMETRIC  PROJECTILES  DUE  TO  MUZZLE  GAS  LOADINGS 


Given  the  physical  characteristics  of  the  projectile  and  propellant 


gas  Mach  number  prior  to  ejection,  V /c..,  the  transverse  linear  an.i 


p'  1’ 


angular  velocities  due  to  both  in-bore  and  muzzle  jet  gasdynamic  load- 
ings may  be  calculated  using  Equations  (31)  and  (32);  however,  the 
momentum  impulse  function  used  must  be  the  sum  of  in-bore  and  muzzle 
jet  impulses,  i.e.: 


P + P, 
o b 


(47) 


The  effect  of  these  transverse  velocities  upon  subsequent  projectile 
motion  will  be  evaluated  in  this  section.  The  angular  deflection,  0 , 

from  the  desired  trajectory  due  to  transverse  linear  velocity  is: 


Qt  = W /V 
t op 


(48) 
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Substituting  Equations  (31)  and  (47): 

0 = (y+1)  p*  S A a — j PT  . (49) 

1 M V Z 1 

P P 

The  effect  of  angular  velocity  is  more  complicated,  requiring  knowledge 
of  projectile  aerodynamics  and  subsequent  oscillatory  motion. 

The  equations  of  motion  of  a statically  stable  missile  have  been 
integrated  by  Murphy  and  Bradley.  They  obtain  the  following  expression 
for  the  aerodynamic  jump,  i.e.,  the  angular  deflection  of  the  trajectory 
which  is  independent  of  the  initial  transverse  velocity: 


0 . = Jzl  Z>'  + J~  l + J * . (50) 

j K 0 £ o A *•  J 

The  first  term  on  the  right  hand  side  is  the  jump  due  to  initial, 
transverse  angular  velocity,  £ ' . The  remaining  terms  represent  the 

o 

contributions  of  initial  yaw  angle,  £ , and  asymmetry  function,  4>.  Murphy  and 
9 0 

Bradley  show  that  for  reasonable  yaw  levels,  the  contribution  of  yaw 
angle  is  negligible,  i.e.,  J~  <<  J~* , Since  the  present  section  con- 
siders only  the  symmetric  missile,  the  last  two  terms  in  Equation  (50) 
will  be  neglected,  thus: 

’ (51) 


where : 


CL 

t~/  _ , 2 ot 

C = kt  c„ 


= Free  flight  lift  and  static  moment 
a coefficients. 


£ = Projectile  shaft  diameter. 

Substitution  of  Equations  (32),  (52)  and  (53)  into  Equation  (51)  yields: 

CL 

Q,  = T2-  (Y+1)  P*  |A  a -A-— — P (54) 

J M M V £ 

a P P 


Of  E<ua™54^“c4™  f ‘°  "UZ:le  £aS  e£feCtS  iS  the  "<*«  “ 


9 - [ 1 + ___  _ ] (y+i)  p*  £ A a _D 

Mx  MV' 

P P 

or  small  angles  of  attack,  the  ratio: 


r 


P^'actu;  inlo^ard^hr161  °f  8r“Vit5'  * c"t*  of  *«««•  “f 

Kith  this  approximation,  the  expression  in  brackets  becomes: 

f 1 * r—  T)  =H  -h  ■ , 


reverse“f  l*nw  ' Jf  the  Center  of  Pressurp  f»r»ard  and 

reverse  flow  coincide,:,  » 4,  the  resultant  jump  of  a statically  stable 

projectile  due  to  muzzlf  effects  is  zero;  houever,  Gallagherll  also 

notes  that  such  a situation  does  not  generally  oc^ur  “nee  the  “Inter 

forward^ light  .USU““y  Cl°S6r  *>  thB  ««  *»  -arse  flop  than  in 

any  fPre?s5-d.  i"  Equation  (55)  may  be  applied  directly  to 

ca  r,,  :;  f-llZed  prol eotrle/gun  configuration.  In  Appendix  B,  sample 
culations  are  presented  demonstrating  the  procedure  to  follow  in 

ZTLlC  reP?n-  aThe  f0r"  °f  Et>UatiOT  (55)  Permits  considimion  of 

fnes  dThe°ie' Ctl1!  d6Slf  A4Crit2ria  t0  minimize  j™P  due  to  muzzle  load- 
ings. The  parameter  p*/Mp  Vp2  may  be  considered  an  indicator  of  the 

ratio  of  propellant  gas  energy  (per  unit  volume)  to  the  projectile 

emcLS  a^C  - RJduction  of  this  parameter  corresponds  lo  a more 

a ?L  rrfer  of  pr°Pellant  energy  the  projectile  in-bore  with 
a resultant  decrease  m muzzle  gas  loadings;  however,  this  also  can  be 


1 +1  G<Zlnl/fnherjr^  E}eme.nts  mch  have  Contributed  to  Dispersion  in 

the  90 /40mm  Projectile,  " BEL  Deport  1013 , U.S.  Army  Ballistic 

Aberdeen  Proving  Ground. , Maryland. , March 


viewed  as  the  standard  "lei.gthen-the-gun-tube-to-decrease-muzzle-blast  " 

CM 

approach  which  has  obvious  physical  limitations.  Increasing  a the 

CL 

a 

static  margin,  also  reduces  jump  due  to  muzzle  effects,  but  may  have 
adverse  effects  on  projectile  interior  ballistics. 


IV.  MUZZLE  BLAST  EFFECTS  ON  THE  FLIGHT  OF  ASYMMETRIC, 
FIN-STABLIZED  PROJECTILES 

In  the  previous  section,  attention  was  given  to  the  deflection  of 
a symmetric  projectile  flying  at  a fixed  angle  of  attack  through  the 
muzzle  gases.  A similar  deflection  would  result  from  an  asymmetrically 
aligned  fin;  however,  in  addition  to  the  jump  caused  by  initial  angular 
velocity,  the  influence  of  aerodynamic  asymmetry  in  free  flight  would 
require  consideration.  Murphy  and  Bradley9  show  that  the  lateral  de- 
flection due  to  aerodynamic  asymmetry  is  dependent  not  only  on  the 
magnitude  and  sense  of  the  asymmetry  ard  the  steady  state  roll  rate, 
but  also  upon  th*.  variation  in  roll  rate  during  spin-up  to  the  steady 
state  condition.  They  demonstrate  that  as  the  launch  roll  rate  de- 
creases from  a steady  state  value  toward  zero,  the  lateral  deflection 
of  the  projectile  increases  significantly.  Their  calculations  only 
treat  initial  roll  rates  which  are  in  the  same  sense  as  the  steady 
state  roll.  Since  the  canted  fins  on  fin-stabilized  projectiles 
launched  from  smooth-bore  o”  low  twist  guns  may  generate  reverse 
(negative)  roll  rates  in  the  muzzle  gases,  the  analysis  of  Murphy  and 
Bradley  will  be  extended  to  consider  this  condition. 

The  jump  due  aerodynamic  asymmetry  is  taken  from  Equation  (50)  as: 


Where 


0.  = J.  *. 

A A 


J a r 

A 2M  lLN 
p e 


CL  CM 
a e 


J e e Te  , 


* Tor  aerodynamic  asymmetry  arising  from  a single  control  surface 
deflected  at  an  angle  r.  with  initial  orientation  angle  <p  . The  ex- 
pression in  brackets  may  be  reduced  to: 

Si  CL 

Me  Cw  " ] ~ Me  1 " S + l)  ' 


Where:  A£  = c.g.  - c.p.  separation  of  asymmetric  fin 
Ac  = c.g.  - c.p.  separation  of  projectile. 

mi  • 4.  ...  x * 


. C A — r J ° 

S>  i-M  t = Af  eA  = 0 which  is  similar  to  the  result  observed  for  the 
lection  due  to  muzzle  blast  effects  on  symmetric  projectiles. 


1Jm  1 rs  ,s2  i^tSj) 

= s+°°  — f f e 1 ds,  ds, 
s o o 1 c 


(60) 


Murphy  and  Bradley  show 

¥ - <j>' 

.*>-■'*  0 


*(s)  M:s+  (- 


Cs 

-)  (e"  -1)  . 


(61) 


Where : 


<t>  • Roll  a.agle  (Figure  7), 

<t>  = Initial,  free  flight  orientation  of  fin  asymmetry  force. 


e = Magnitude  of  fin  cant. 


= Initial  roll  rate  (after  muzzle  blast), 
C, 


'A,  6 

j « 6 


Ca  + k 2 c 

p a D 


- Steady  state  roll  rate. 


, I 
2 x 


» ‘7? 

P 


C = - ('2~'^~')  (^a2  C£  + cq)  * Roll  damping  coefficient, 
P a P 


Roll  moment  coefficient  due  to  roll. 


C = Roll  moment  coefficient  due  to  cant. 


6 = Cant  angle. 


To  obtain  a generally  applicable  solution  to  the  integral  in 
Equation  (60),  the  following  non-dimensional  parameters  are  defined: 


st  = Cs, 

a! 


A = 


T » 

a1 


B = -2L 


♦'  . 


Thus,  Equation  (61)  becomes 
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(62) 


<f>  (^0  = A [ st  - (B-l)  fe 


and,  substituting  into  Equation  (60) : 


lim 


® =(s. 


) 


Cs. 


f 

0 


I 

0 


'2  i<Ks.  /C) 

6 tl 


Murphy  and  Bradley  show  that  for  constant  rolling  motion. 


d s d s 
1 z2 


• (63) 


$ = 


(64) 


and  use  this  value  as  a standard  against  which  to  index  the  magnitudes 
of  jump  with  varying  spin: 


Thus : 


t » 


t c 

lim  , s t„ 

♦t - A(st  ”)t!  { d » a s 

t 1 r2 


(65) 


(66) 


The  double  integral  in  Equation  (66)  is  evaluated  numerically 
using  a Runge-Kutta  technique  with  various  values  of  A and  B.  The 
resulting  values  of  the  magnitude  of  $ are  plotted  versus  4>^/C  for  various 

» Figure  16.  The  abscissa,  is  a measure  of  the  angle 

turned  through  as  the  rol]  rate  approaches  the  steady  state  value.  The 
behavior  of  |^^.|  indicates  that  as  the  number  of  revolutions  required  to 

achieve  steady  state  spin  increases  (i,e,,  increases),  the  jump 

due  to  asymmetry  also  increases.  The  current  technique  is  seen  to  pre- 
dict values  of  | | which  are  identical  to  those  calculated  by  Murphy  and 

Bradley  for  = 0,0.1  and  1.0.  The  curves  clearly  show  that  as  the 

initial  roll  rate  decreases  from  the  steady  state  value,  = 1.0, 

the  jump  due  to  asymmetry  increases  significantly. 

The  orientation  of  0^  is  plotted  in  Figure  17.  Since  the  jump  is 

determined  by  the  product  of  JA  and  $,  Equation  (58),  the  direction  of 

jump  is  simply  the  sum  of  the  two  arguments.  The  Y and  Z components 
resulting  from  the  evaluation  of  Equation  (63)  as  the  limit  is  approached 

are  plotted  in  Figure  18  for  <J//C=31.  This  plot  reinforces  the  conclu- 
sions drawn  from  the  previous  data  regarding  the  behavior  of  asymmetric 
jump  as  the  initial  roll  rate  decreases.  In  Appendix  C,  the  results 
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of  this  section  will  be  applied  to  the  calculation  of  asymmetry- 
induced  jump  of  XM-645  flechette. 

V.  CONCLUSIONS 

An  analysis  is  developed  permitting  the  computation  of  the  magni- 
tude of  gasdynamic  loadings  on  fin-stabilized  projectiles  during  launch. 
Both  in-bore  loadings  due  to  the  unsteady  expansion  propagating  upstream 
and  external  loadings  generated  in  transit  of  the  muzzle  blast  are 
considered.  The  in-bore  flow  is  modeled  as  a one-dimensional,  unsteady 
expansion  fan,  while  the  external  flow  is  approximated  as  a quasi-steady, 
underexpanded  jet.  A method  of  characteristics  code  developed  by  NASA 
is  used  to  calculate  the  external  flow  properties. 

The  transverse  loadings  on  the  projectile  are  assumed  to  be  dom- 
inantly generated  on  fin  surfaces . As  such,  they  are  calculated  using 
two-dimensional,  thin  airfoil  theory.  The  nature  of  the  assumptions 
made  to  simplify  the  analysis  results  in  a solution  which  is  generally 
applicable  to  any  projectile  fin  configuration  provided  the  sabot  design 
is  a center  or  forward  puller.  Weapon  characteristics  enter  the  solu- 
tion in  a straightforward  manner  permitting  direct  computation  for  any 
caliber,  launch  velocity,  or  exit  pressure.  Additionally,  the  assumptions 
produce  an  upper  bound  on  the  magnitude  of  muzzle  blast  loadings.  This 
permits  the  simplistic  approach  contained  in  this  report  to  be  used  to 
estimate  projectile  dispersion  due  to  muzzle  blast  loads.  Comparison 
of  this  prediction  with  observed  values  of  dispersion  will  demonstrate 
the  relative  importance  of  gasdynamic  effects. 

The  analysis  of  the  in-bore  flow  is  highly  idealized;  however,  it 
provides  an  estimate  of  the  magnitude,  duration,  and  relative  importance 
of  loadings  seen  by  the  projectile  in  this  region.  The  momentum 
transferred  to  the  projectile  is  shown  to  be  a function  of  propellant 
gas  properties  prior  to  separation  of  the  obturator  and  the  standoff 
distance  between  the  obturator  and  the  fins.  For  high  velocity  guns, 
in-bore  loadings  may  be  neglected;  however,  for  low  velocity  guns, 
such  as  mortars,  they  must  be  considered. 

The  method  of  characteristics  calculation  of  the  external  flow 
used  in  this  report  is  shown  to  produce  momentum  impulses  significantly 
greater  than  those  previously  computed  using  source  flow  models. 
Comparison  of  the  momentum  impulse  generated  under  quasi-steady  and 
fully  steady  approximations  of  the  muzzle  jet  flow  demonstrate  that 
adequate  results  are  obtained  with  the  simpler  steady  analysis. 

Further,  sample  calculations  for  an  XM-645  flechette  are  presented 
showing  transverse  momentum  imparted  in  the  muzzle  blast  to  be  an 
insignificant  contribution  to  the  total  observed  dispersion  of  the  round. 
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8 = FIN  CANT  ANGLE 

r0=  DISTANCE  FROM  PROJECTILE  AXIS  TO  FIN  C.P. 


A.  DIFFERENTIAL  FIN  CANT 
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a = PROJECTILE  ANGLE  OF  ATTACK 
A = DISTANCE  FROM  PROJECTILE  C.G.  TO  FIN  C.P. 

B.  PROJECTILE  AT  ANGLE  ATTACK 

Figure  13.  Projectile  Launch  Configurations 
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ANGLE  OF  ATTACK 


Consider  a multi-finned  projectile  at  angle  of  attack  relative  to  a 
reverse  flow,  V,  Figure  Al.  The  attitude  of  fin  A (in  the  plane  per- 
pendicular to  the  plane  formed  by  the  projectile  axis  and  the  velocity 
vector)  is  identical  to  that  of  the  projectile  axis  and  is  determined 
from 

v * n = cos  3 = -sin  5.  (Al) 

Where  v - Unit  vector  in  direction  of  V, 


n - Unit  vector  perpendicular  to  surface  of  fin. 


Since 


n = i * 

(A2) 

v = cos  3 i - sin  3 j . 

(A3) 

for  small  3, 


1 

;; 

i 


V « 1 


a j 


(A4) 


and 


v*n  = - a. 


(A5) 


It  is  of  interest  to  consider  a fin  at  arbitrary  orientation  on  the 
projectile  body.  Figure  A2.  The  angle  between  the  fin  and  the  axis  is 

2irk 


(A6) 


Where 


From  the  diagram  f 


4>0  - Angle  between  j_  axis  and  first  fin, 
n - number  of  fins, 
k - fin  index  number. 


sin  <j>  j - cos  <J>  k 


(A7) 


The  lift  on  this  fin  may  be  obtained  from  Equation  (8) : 


1 2 

£ = - C£  3 sin  $ ^ PV  A . 

a 


0*9) 


For  small  3,  the  components  of  the  force  in  the  j,  k,  plane  are, 


£ = - C£  3 i pV2  A sin2  <f, 

?v  a 

1 2 

£ = C,  2 r pV  A sm  <f  cos  <f. 

4 

Yh  a 


(A10) 

(All) 


The  total  vertical  and  horizontal  components  of  force  exerted  on  the. 
projectile  due  to  the  combined  loadings  on  the  fin  array  may  be  obtained 
from  a summation  of  Equations  (A10)  and  (All)  : 

n 1 ,,2  , . 2 ,,  x 27Tk., 

L = 2 ( -C.  a pV  A sm  (<f  + — — ■)] 

v k-1  a * 


1 2 . . 2 2irk. 

= - C 3 j PV  A 1 sin  C*o  + IP  ' 


“2 
a 

, 1 


k=l 
n 


Lh  = CjC  a 7 pv^  A s sin  + 5 cos  »o  + PP1  ' (A13) 

a k*l 

Using  trigonometric  identities,  the  summation  on  the  right  hand  side 
of  Equation  (A12)  may  be  expanded  to 


(A12) 


2irk. 


n 

2 

k=l 


n 


. 2 . 2Ttk,  _ ? .2  , U — 2 2 rrk  . „„2  , ? _,_2  2Trk 


2 sin^  (<fQ  + — ) = sin  <fQ  2 cos  — 

K — 1 


+ cos  (j)  Z sin 
0 k-1 


n 


, Z ■ 47rk 

+ sm  if  cos  if  2 sm  — 

0 k-1  n 


(A14) 


Jolley12  gives  the  following  sums  for  the  above  series: 


12.  L.  B.  W.  Jolley,  Summation  of  Series , Dover  Publications,  New  York, 
N.Y.,  1961.  
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I sin2  kQ  = 2.  - £PJ-L^e^siiniin8_ 

^ 2 2 w»1T1  0 


n 

2 sin  k0 
k=l 


sin  (n+1)]  0 sin  i n3 


G 


sm  7 


2tt 

For  0 = — , Equation  (A15)  shows 


n 


_ . 2 , 2v  n * 

2 Sin  k — = rr  , 

k=l  n 2 


and  since: 


" . 2 2TTk  " 2 2irk  " . 2 2TTk  2 

2 sin  + 2 cos  — r—  = 2 (sm  — — + cos 


k=l 


n 


k--l 


n 


k=] 


n 


= n. 


Equation  (A17)  and  (A18)  show: 


n 

2 cos 
k=l 


2 2 it 


n 

2 


4tt 


For  0 = — , Equation  (A16)  shows 


4ir 

2 sin  k — 
k=l  n 


= 0. 


0 

*This  is  valid  for  n >_  3.  For  n = 2,  sin  - sin  0 = 

an  indeterminate  form  on  the  right  hand  side  of  Equatio 

n = 1,  the  need  for  a summation  vanishes  and  L ~ sin^ 

’ u 

cos  <j>  . 
o 
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2trk' 

n 


0 yielding 
n (A15) . F 

*0-  ' si 


p,  map  pmihl.jp. 


Substituting  Equations  (A17) , (A19),  and  (A20)  into  Equation  (A14) : 


v e.-  2 , 27rk.  n . 2 n 2 , 

Z sin  (4>  + •)  = ■=•  sin  $ + cos  6 

^ 1 o n ' 2 To  2 ' o 


n 

2 ‘ 


(A21) 


The  resulting  expression  for  the  vertical  component  of  lift,  Equation 
(A12),  becomes 

-v  (A22) 


t r*  v 1 2 . n 

K * - C£  3 7 PV  A J 

a 


Using  a similar  procedure,  it  may  be  demonstrated  that  the  summation  in 
Equation  (A13)  is  equal  to  zero;  thus 


Lh  = 0 . 


(A23) 


Equations  (A22.)  and  (A23)  indicate  that  the  lift  on  a multiple  fin 
assembly  (n  >_  3)  is  independent  of  roll  angle,  <j>^,  and  acts  in  the  plane 

formed  by  the  flow  velocity  vector  and  projectile  axis.  The  loading  on 
an  n-finned  assembly  is  equivalent  to  that  on  a two  fin  assembly  with 
its  fins  oriented  perpendicular  to  this  plane  and  with  a total  fin  area 


of  ^ A. 


; i 


For  the  XM-645  flechette,  some  of  the  parameters  necessary  to  complete 
this  calculation  are  given  in  Table  I.  The  values  of  p*  and  require 

information  on  the  emptying  of  the  gun  tube.  The  exit  properties  of 
the  5.77mm  smoothbore  gun  firing  this  flechette  at  1462  m/s  were 
calculated  by  A.  Celmins  of  BRL  using  a technique  described  in  Ref- 
erence 13. 


The  pressure  and  Mach  number  at  the  muzzle  obtained  from  these 
calculations  are  shown  in  Figure  Bl.  The  behavior  of  the  pressure  is 
typical;  however,  the  variation  of  Mach  number  with  time  is  somewhat  un- 
usual. According  to  Celmins,  the  Mach  number  behind  the  projectile  is 
supersonic  prior  to  shot  ejection;  further,  the  muzzle  Mach  number  in- 
creases from  1.23  as  the  fins  pass  the  muzzle  to  1.36  as  they  exit  the 
supersonic  jet  core.  During  the  same  period,  the  muzzle  pressure  de- 
cays from  382  atmospheres  to  363  atmospheres.  The  impact  of  temporal 
variations  in  muzzle  conditions  on  the  analysis  presented  in  the  body 
of  the  report  will  be  addressed  in  these  calculations. 

The  quasi-steady  approximation  of  Oswatitsch1  postulates  that  jet 
properties  may  be  computed  throughout  the  flow  field  based  upon  the 
muzzle  exit  conditions  at  the  time  of  interest.  Using  this  approach, 
the  jet  properties  at  any  instant  during  gun  tube  emptying,  are  defined 
by  interpolating  the  centerline  Mach  number  distributions  computed  by 
the  method  of  characteristics.  Figure  6B,  to  match  the  muzzle  exit  con- 
ditions, Figure  Bl.  As  the  fins  traverse  this  time  varying  flow,  the 
lift  function  may  be  calculated  based  on  the  instantaneous  flow  proper- 
ties and  Equation  (12): 


To  evaluate  the  momentum  function,  it  is  necessary  to  include  the  vary- 
ing p*  in  the  integration  in  Equation  (23) : 


★ • 

Where:  p^  = critical  pressure  just  prior  to  shot  ejection 

★ 

p = critical  pressure  varying  with  time. 

The  time -dependent  momentum  function  becomes: 

P.  , = /*  K L dx  . (B2) 

t.d.  o p* 

1 

Both  the  lift  and  momentum  functions  have  been  numerically  evaluated 
and  are  shown  in  Figure  B2.  The  asymptotic  value  of  ?f . is  0.287. 

Substituting  into  Equation  (55)  : 

0 = - 0.078  d .mils (for  a in  degrees).  (B3) 

This  deflection  has  a magnitude  of  0.078  mils  for  each  degree  of  pro- 
jectile angle  of  attack  at  launch.  The  deflection  is  directed  opposite 
to  the  angle  of  attack  orientation. 

As  developed  in  the  body  of  the  report,  the  analysis  of  muzzle 
blast  effects  on  the  projectile  trajectory  assumes  steady  muzzle 
properties  throughout  projectile  residence  in  the  supersonic  core  of 
the  jet.  The  value  of  deflection  predicted  by  the  steady  model  will 
be  compared  with  the  time-dependent  result,  Equation (B3)  . Assuming 
muzzle  exit  conditions  remain  steady  at  the  values  seen  by  the  fins 
as  they  pass  the  muzzle  (Mg  = 1.22,  pg  = 382  p^) , the  momentum  impulse 

during  jet  residency  is  0.295,  Figure  12.  Substituting  into  Equa- 
tion (55)  : 

0 = -0.091  3 mils  . (B4) 

The  steady  deflection  is  seventeen  percent  greater  than  the  value 
given  in  Equation  (B3) . This  overestimation  would  be  anticipated  since 
the  unsteady  calculation  reflects  the  decay  of  momentum  flux  from  the 
muzzle  as  the  gun  tube  empties,  while  the  steady  results  do  not.  The 
fact  that  the  steady  approach  results  in  a slight  overestimate  is 
valuable  in  that  it  permits  a straightforward  computation  of  the  upper 
bounds  of  muzzle  blast  induced  jump.  If  this  jump  is  significantly 
less  than  the  observed  total  dispersion  of  the  round  in  question,  muzzle 
gas  effects  may  be  discounted  as  an  error  source,  eliminating  the  need  to 
calculate  the  fully  time  dependent  flow. 

Schmidt  and  Shear5  measure  launch  angle  of  attack  of  the  XM-645 
to  be  in  the  range: 


. IL 


Substitution  of  this  spread  of  launch  attitude  into  Equation  (B3)  results 
in  a jump  range: 


0.055  mil-  < 0 < 0.016  mil 


Comparison  of  this  magnitude  of  jump  with  recorded  dispersion  of  the 
XM-645  and  related  flechettes^  (on  the  order  of  one  mil  or  greater) 
establishes  that  muzzle  blast  does  not  contribute  significantly  to  the 
total  dispersion  of  the  round.  Other  sources  of  transverse  jump,  such 
as  in-flight  asymmetry  and  in-bore  vibration,  must  be  examined. 


J.  F.  Thompson^  'Evolution  of  the  Sabot  Design  for  the  Serial 
Fleohette  Rifle 3"  U.S.  Army  Small  Arms  Systems  Agency 3 Aberdeen 
Proving  Ground , Md. , May  1972. 
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Muzzle  Exit  Properties  of  5.77mm  Smoothbore  Firing 
XM-645  Flechette 


Figure  B.l 
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APPENDIX  C:  SAMPLE  CALCULATION  OF  MUZZLE  BLAST  EFFECTS  ON  AN 

ASYMMETRIC  XM-645  FLECHETTE 

The  total  jump  of  an  asymmetric  projectile  may  be  expressed  as  a 
vector  sum  of  the  jumps  due  to  asymmetry,  initial  transverse  angular 
velocity,  and  initial  transverse  linear  velocity: 


6 = JA  0 + V f 


To  isolate  the  jump  due  to  aerodynamic  asymmetry,  the  projectile  con- 
figuration shown  in  Figure  Cl  will  be  considered.  The  round  is  assumed 
to  be  launched  with 


“i- 


= 5 0 


The  asymmetry  is  postulated  as  two  opposing  fins,  inclined  with  an 
angle  e (in  the  same  sense  for  both  fins)  with  respect  to  their  normal 
orientation.  The  differential  fin  cant  angle,  2o,  between  opposing  fins 
is  unaltered;  thus,  the  steady  state  roll  rate,  is  maintained. 

Additionally,  the  orientation  of  the  fins  at  launch  is  in  the  X , Y 

6 G 

plane.  Figure  Cl. 

The  transverse  angular  and  linear  velocities  imparted  by  the 
muzzle  jet  may  be  calculated  using  the  results  of  Appendix  B.  Since 
the  present  technique  only  considers  fin  loadings,  the  impulses  im- 
parted to  the  asymmetric  fin  pair  at  an  angle,  e,  are  identical  to 
those  experienced  by  a symmetric  (four-finned),  XM-645  projectile 
at  an  equal  angle  of  attack,  e.  The  jump  due  to  the  transverse 
loadings  may  be  taken  directly  from  Equation  (B4)  as 


Jj  ?/  «.  JL 

\ Vr 


-0.091  e e1  ^e+1T) 


where  e = Magnitude  of  asymmetric  fin  cant, 

4>e  = initial  orientation  of  asymmetric 
force  in  free  flight, 

<f>  = orientation  of  asymmetric  force  in 

the  muzzle  jet. 

Calculation  of  the  jump  due  to  aerodynamic  asymmetry,  UA$,  is 

more  complicated  since  it  requires  evaluation  of  the  reverse  spin 
imparted  in  the  muzzle  jet  and  estimation  of  certain  properties  of 


65 


the  XM-645  not  given  in  Table  I.  The  terms  in  the  expression  for 
asymmetric  jump  may  be  given  for  this  geometrical  configuration  as 


r £ 

2M  1 N 


Cl  Si 

-2—^  ] e el+e  , 


* = \/^  , (C4) 

= ®t  C+l/C.+'/O  » (Defined  in 

Figures  16  and  17)  (C5) 

C =*5r1  (ka"2  C£  + V * (C6) 

P P 

The  normal  force  coefficient,  CN  , is  defined  as 

e 

L 

CN  = 1 (C7) 

e ipV^  Se 
2 P 

where  = Lift  on  projectile  due  to  the  asymmetric  fin  pair. 

The  lift  may  also  be  evaluated  in  terms  of  the  lift  force  acting  on  the 
asymmetric  fin  pair,  expressible  as 

L e = 7 p V 2 (2A)  C£  e . (C8) 

Substitution  of  Equation  (C8)  into  Equation  (C7)  yields 


C = C 
N S jC 
e a 


For  the  XM-645, 


Mm  = 4.36; 


Thus,  according  to  thin  airfoil  theory  the  lift  coefficient  is 


C.  = 4 (M2  -I)'3*  = 0.942  , 


(CIO) 


- -Ill 


Wl 


I 


, . » ^ J 


permitting  the  evaluation  of 


,-l. 


CN  =7.86  (rad  ) 
e 


Ceil) 


For  slender  projectiles  at  small  angle  of  attack: 


"M 


JN 


€ 

£ ’ 


(C12) 


where 


A = c.p.  - c.g.  separation  of  asymmetric  fin 
= 9.9  £ (assuming  c.p.  at  fin  centroid). 

Fquations  (Cll)  and  (C12)  permit  the  calculation  of 


-2  i4>f 

J = -0.11  X 10  ee  mils(for  e in  degrees).  (C13) 

i\ 


The  determination  of  the  value  of  $ requires  values  for  the  initial 
(subsequent  to  penetration  of  the  muzzle  gases)  and  steady  state  roll 
rates,  4/  and  4>'  , respectively.  Murphy  and  Bradley^  give  the  following 


relation: 


= - 


'l6  5 


C*  tka2cD 
P 


C14) 


which  may  be  rearranged  and  substituted  into  Equation  (C5) , the  expression 
for  C : 

C„ 


c = (_4 — ■) 


2 M \ 2 
P ka 


(Cl  5) 


The  roll  moment  coefficient  due  to  fin  cant  is  defined  as 


M 


(C16) 


i P V*  S£S 

Z p 


Again,  the  roll  moment  may  be  expressed  as  a function  of  the  lift 
on  the  differentially  canted  fins. 


Mx  - \ 0 vp2  '4A>  T CX  S 
r a 
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ESSJ 
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(CT7) 


M 


: 


where  the  moment  arm  from  the  projectile  axis  of  symmetry  to  the  fin 
center  of  pressure  is  assumed  to  be  D/4.  Substitution  of  Equation  (C17) 
into  Equation  (C16)  yields 


AD 


r - n r 
\ ' 51 


.-1, 


t - ST  ’ 1U89  (rad  } 


(CIS) 


Equation  (C18)  and  Table  I are  used  to  evaluate 


and 


C = 5. 89  X i(J  (per  caliber), 
<j/  /C  = 23.8  radians. 


(C19) 

(C20) 


Since  the  launch  roll  rale  is  taken  to  be  zero,  the  roll  rate 
upon  penetration  of  the  muzzle  gases  may  be  computed  using  Equation  (28), 
which  shows  that 


nr 


*/  = 2.  <y+i)  p A 6 P 

Vo  I 2 ''  } v ( 

x V 

P 


The  previous  Appendix  demonstrated  that  reasonable  accuracy  is  obtained 
if  variations  in  muzzle  conditions  with  time  are  neglected:  thus. 


assuming  M = 1.22  and  p = 382  atmospheres: 

G ® 


d/  = -5.65  X 10  rad/cal, 
o 


(C21) 


and 


/*•'  = -0.0403 


(C22) 


From  the  values  of  <t>^/C  and  <)/ /<f>^  in  Equations  (C20)  and  (C22) , 


$ may  be  determined  from  Figures  (16)  and  (17)  to  be 

t ' fdl0! 

* * 7.13  e11  J (radians), 


and 


* - 509  e1*41  > 
Equations  (C13)  and  (C24)  yield 


(C23) 

(C24) 
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(C25) 


* 


J,  4 = 0.56  e e 


(*E  * 


TT  + 


41°) 


mils  (for  e in  degrees) . 


It  is  now  possible  to  evaluate  the  total  jump,  Equation  (Cl),  by 
summation  of  the  jump  due  to  transverse  velocities.  Equation  (C2) , and 
the  jump  due  to  asymmetries,  Equation  (C25) . This  summation  is  depicted 
graphically  in  Figure  Cl,  and  the  result  given  below 


i (138°) 

6 = 0.49  e e mils  (for  e in  degrees).  (C26) 


An  interesting  comparison  may  be  made  between  this  value  of  jump 
(including  muzzle  blast  effects)  and  a value  computed  under  the  assumption 
that  the  launch  properties  of  the  projectile  are  unaltered  in  transit 
uf  the  muzzle  gases. 


For  the  launch  conditions  considered  in  this  calculation. 
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the  total  jump  (in  the  absence  of  muzzle  gas  effects)  is 


0 = 0A  = JA  * 


(C27) 


JA  is  given  in  Equation  (C13)  and  $ may  readily  be  determined  knowing 

t /C  = 23.8, 
and 


I J 

' I 


* I 


0.499  e e ^ ^ mils  (for  e in  degrees)  (C30) 

Comparison  of  Equations  (C2S)  and  (C30)  shows  that  the  effect  of 
reverse  spin  imparted  in  transit  of  the  muzzle  jet  is  to  increase  the 
magnitude  of  jump  due  to  asymmetry  by  12  percent  over  the  value  for 
zero  spin.  Equation  (C26)  shows  that  this  increase  is  nearly  exactly 
cancelled  when  summed  with  the  jump  due  to  transverse  velocities  genera- 
ed  in  the  muzzle  jet.  Not  only  are  the  magnitudes  of  jump  given  in 
Equations  (C26)  and  (C30)  almost  identical,  but  the  orientation  of  jump 
is  within  one  degree  for  the  two  cases.  Thus,  for  the  XM-645 ' rw.zzle 
blast  loadings  do  not  significantly  alter  the  jump  experienced  in 
launching  an  asyrmetric  projectile.  This  conclusion  is  of  course 
qualified  since  only  one  projectile  and  set  of  launch  conditions  were 
considered;  however,  the  result  seems  to  justify  a more  general, 
mathematical  analysis  of  this  effect. 
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PHYSICAL  PROPERTIES 


DYNAMIC  PROPERTIES 


0.014  rad/ cal 


TABLE  I:  PROPERTIES  OF  XM-645  FLECHETTE 


LIST  OF  SYMBOLS 


planform  area  of  single  fin 

local  sound  speed 

roll  damping  coefficient 


roll  moment  coefficient  due  to  roll 


roll  moment  coefficient  due  to  fin  cant 


CL’  CL 


projectile  lift  coefficient  and  derivative  with  respect  to  a 


C£  ’ C£ 


fin  lift  coefficient  and  derivative  with  respect  to  a 


projectile  static  moment  coefficient 


moment  coefficient  due  to  fin  deflection 


normal  force  coefficient  due  to  fin  deflection 


distance  from  obturator  to  fin  c.p. 
diameter  of  gun  bore 


jaj  Ji‘  Ji' 


longitudinal  and  transverse  moments  of  inertia 


jump  coefficients 


V kt 


axial  and  transverse  radius  of  gyration 


shaft  diameter  of  projectile 
projectile  and  fin  lift  forces 


dimensionless  fin  lift  function 


Mach  number  or  moment 


mass  of  the  projectile 


number  of  fins 


LIST  OF  SYMBOLS  (Continued) 
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P 

P,  P 
P 


P 

t 

s 


5 
t 

V 

V 

p 

V 

r 

w 

X,Y,Z 

a,  g 

V 

6 

A,  Af 

€ 

l 

0 

6t 

0. 
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pressure 

momentum  transfer  to  projectile  and  single  fin,  respectively 

dimensionless  momentum  function 

muzzle  jet  and  in-bore  momentum  impulse 

total  momentum  impulse,  P + p 

o b 

distance  along  trajectory  in  calibers  (s  = x /£) 
reference  area  of  projectile  (S  = 7t£2/4) 
time  (t-0  corresponds  to  obturator  passing  muzzle) 
local  flow  velocity  in  muzzle  flow 
projectile  velocity 

relative  velocity  V - V 

P 

transverse  velocity  in  Z direction 

e 

coordinates 

angles  of  attack  and  sideslip  in  non-rolling  coordinate 
system 

ratio  of  specific  heats 
differential  fin  cant  angle 

C-P*  - c.g.  separation  in  reverse  and  forward  flow,  respectively 

magnitude  of  asymmetric  fin  cant  angle 

complex  angle  of  yaw  B + i 3 

angular  deflection  of  projectile  from  boreline 

angular  deflection  of  projectile  due  to  transverse 
velocity  in  muzzle  blast 

aerodynamic  jump:  angular  deflection  from  particle 

trajectory  (gravity  and  drag  determined)  due  to 
aerodynamic  forces 
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denotes  earth-fixed  coordinates 

denotes  conditions  evaluated  at  the  relative  velocity, 

denotes  conditions  in  gun  bore,  at  the  muzzle,  just 
prior  to  obturator  separation 

denotes  projectile  properties  immediately  after 
penetration  of  the  muzzle  blast 


denotes  ambient  or  steady  state  conditions 


denotes  dimensionless  quantities 
denotes  time  derivative 


( )'  denotes  derivative  with  respect  to  s 


I 
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